,rV ! ; V/i' 

VT ^ > _ v jV' ' y J> ' u' -*-s 

v '< ‘ * V '/• 

■ K.; A " ;.W ,yV ( / 

" vV v .^ 7 ,y' 

,v *0, H ^ 

■ ^ 

\ •.'■£/" ‘A- : \ ^ 

: • ' v -a "A..- 

-•■ WW r ‘V 

V V , - f /A\V 

■ "M Xr V ' 

T . ,-N 

. x ^ >. : 


;W v via 


V'W 

-. > 3=^ 




X- 661.- 7 2-447 


s\i.* ' 

• 'v-_ */ 






M VAV' ,.%> 


'Ml.: 


r V ■ Ar 1 

‘TP'S- \t '■* 


vA 'L 


.-"V'' 


- ->- -C 


. i, -0 /" 

‘ ^ Z' 

7 C 7 S • 


- 4 >rj. 7 : >' v V> 7 r Vy '7 T< v 

RESULTS ON THE ENERGY ' 

DEPENDENCE OF COSMIC RAY Vi 
: ,V /“UADftF COMPOSITION 


VA^mwb'oI’coskc^^^ 

V, - (COMPOSITION t *i I CH «R0B 


W ‘ '^>1 

y n ■£ v .,ry 


m ' ;yy 

V - ' l> 

) J > '/• -VU 


UJ? . COSMIC 

COMPOSITION (NASA 1 Jt* CHARGE 

-j uw (nasa) p Ac $ 3 .qo 

i ^ CSCL 03 B 


' A . 

TJ .. 5 1 ‘ 


. - Ai f - : f 7 yA'*> 7 £~~r 

,v ■/-: Aiv w > ' v • < 

A" ■ ■ • A"- 5 * f* > . ' Sv jv V *- M ' 1 


Vv K; BALA SUBRAHM ANY AN 

•• VA. 


/ v -SI' j N - ' !, 


UJe r „ Unclas — 
_ G3/29 06241 h v 

7 v-ii 

H MANY AN r vi 


> ■ 


^7 ' > 


‘i^Tr ’->L ■ . ( : {. - f? y 

In ‘ ^ v . ' — 4 v . > M I. " n‘ ‘.n 

" "ZM* 


A ; A ^" : V V 

n< ' --'•: ‘ 


./ < 

;Vr : / 'A 


• ' : vvA , A A'* V: 

ASV:- 


' ' ^ X- ' AW A — x 
w ■ C ' -A- v- 

-A A ,/'7 A' - ( ; N 

< M/ ^ f l y A-^’v / 

' r- 7 N ~ >>7 ^ r- 

' •• y ' 

^ / A- X \ l N-** \ 


^ r 

< L . . >/ 

A aT v 1 r ' ' V 
j '~ M v -_ aW _ 


' / 


-•> ’ M 

t '-7 


k -M 

7 , ' / < 7 


' ; ‘ V ^ v- 
A-n ■ ■ -M 


,y\ ' • ■•-; V . * 

:-'V , /. , . V 

' W " T- • 4 


Vr;V--iV" 


- : N" : • ' vvi— 

> (N- ~ 1 ' ■ *' ‘\ — \ ' J •■* ^ 1 a y J- 

/ j " • v ^ .;-/Vv v„ v • • y ' • - - 

*■.-. • ' •. ... < ' / V-^. ■ .. 

, ^ ^ y. • V /- . ^ . .^r \ j * 


■ s '-n Vv , , / . A; . 

■ • -DECEMBER 1972 

( Revised mayai973 

I , A-, - \ J ' ^ . ■-- | ‘ 


r 

* v 

< ■' . 

M J- f 



s - k '> 

■r " 1 * 

V. A 

✓ t, w 

'•' V ' 

y ; 

V.A 

.;■> X J 

r c. 

' "■• . 

X • , . A 


•a7 

" J - o- 

AV ;- ^ 


':<t' 

77/^ 

6>y 


. v- ,v .;- 
v./- V ' 


v 7 ::-.-k. 


,; -^r. j ■ V-; 

■//VW V> --VA:; 

-r-rn i^L *■' , A' v \ 7 


/ ^ v 
/ ^ vv 


A.‘^ . , - A 

, /v v V* i 

.; 7 • • ' 

r- V 


;^r : 7\ 


,v ' y \ ' vp 

^ ^ ^ ‘ ’ VV ^ ■ /S? 


v^vvJI 


:> ’ ■ v; 

•A V 


Vptvfl 

- v . /5s? y w • w, r- ^ 

/'- MSS. %‘%: %'a 


r ^ 

GSFC 




f .si' A Y/§ N %%' %> .% 

v ' ■ v i£ A\: •; N w 

^ \ --/ j V§- V - ,vA A 


A VA-A, AV f :- A. A- V- u V^%>V :f - Vy 
' r W ■ A ■'•■' V < -v .-i ; ^ ^>Mi L n L UW 

GODDARD SPACE FLIGHT CENTlf F^ 4 

J •*■ ' v; • V , vA-‘VA GREENBELT, MARYLAND A’ Vf - 

A ~N -K: .‘ -AA • • VAV.; A ^:VVj;, >v ' V A 

A ' - 7 ■> "v ■ '•->'•:./* ';-• .s .NV, ./j. aA ; }'-/' v '<v ._ .. '.a VJ 
v .• ' -. Av . •A <- M . ■: C ,a, •/ c-v ' • i r V ‘A:A V , 

W‘ '' — "' Si '7 • />•:. 'Y • - - itv. • V ^ ^ -V " > • ’< 


bUllUAKU i 

. ^ V.: 77 ;V V" 6REB 
AA- avAv/a 


'/ y:: 7 ^ - • ■ 

w-. .' . c M: \ 

. • ■' 7 " 7 '" - ' '7 • 

C > ' ■ v-j ^ tn f 


in i ’ util itn.^ " 4 Y'- 

rtANfD rT 

■ . v^. '■ /yJ 

<" <\ Mt ^ " - - * 1 7 

- 7 -'-:(r;;; V Av 7 : 7 ^v: 
: / V Yv v 77 -v..r-^ -^ N 
V' r 7 Txy y 



RESULTS ON THE ENERGY DEPENDENCE OF 
COSMIC RAY CHARGE COMPOSITION 


V. K. Balasubrahmanyan and J. F. Qrmes 
Goddard Space Flight Center, Greenbelt, Maryland 20771 

ABSTRACT 

Measurements using a balloon borne ionization spectrometer on the 
differential energy spectra of the heavy nuclei of the galactic cosmic 
radiation are reported. The results include more data and improved 
charge and energy analysis over our previously reported results. The 
spectra of individual elements up to oxygen and groups of nuclei up 
through iron have been measured up to almost 100 GeV/nucleon. The energy 
spectrum of the secondary nuclei , B+N, is steeper than that of the primary 
nuclei, C+0, by y = 0.21 + .09 in agreement with Smith et al., 1973. The 
spectral shapes found by us are reasonably well represented by single 
power laws between 2 and 60 GeV/nucleon. Our data are consistent with the 
decrease in the secondary to primary ratio found by Julius son et al. (1972) 
above 20 GeV/nucleon, but our data show no evidence for any sudden change 
in this ratio within counting statistics. The most dramatic finding is 
that the spectrum of the iron nuclei is flatter than that of the carbon 
and oxygen nuclei by 0.57 +0.14 of a power. The 10 <. z < 14 group spectrum 
is consistent with that of C+0 within errors. The experimental techniques 
for charge and energy determination are presented. Charge resolution is 
unique through oxygen, and at iron the error is + one charge, so the 
charge groups are clearly separable. Energy is measured with a - 30% 
between 2 . 5 and 25 GeV/nucleon . Corrections ‘due to nuclear disintegration 
and losses of energy out the bottom of the spectrometer are discussed. 



I. INTKDDUCTION 


Until recently the detailed study of the energy distribution of galactic 
cosmic rays above 10 GeV/nucleon was confined to integral measurements using 
threshold devices (von Rosenvinge , 1970, Webber et al., 1973) . Current ex- 
perimental programs of the Berkeley group (Smith et al., 1973), university 
of Chicago (Juliusson et al. , 1972) and Goddard Space Flight Center (Ormes 
and Balasubrahmanyan, 1970) are attempting to extend the knowledge of charge 
composition to approximately 100 GeV/nucleon. For the first time differential 
measurements of the energy spectra of individual charges and charge groups 
above 10-GeV/nucleon are being reported (Ormes et al., 1971, and Smith et al., 
1973) . The new results in this energy region point towards very significant 
developments in our understanding of the life history of cosmic rays. Below 
a few GeV roost of the experimental results were consistent with an energy 
independent composition. After a thorough analysis in which an improved 
charge resolution was coupled with detailed trajectory and energy analysis, 
we conclude that the cosmic ray composition varies with energy above 1 GeV/ 
nucleon. This will have important implications for cosmic ray astrophysics. 

It is only by studying all of the relevant indicators, the seoondary/priirary 
ratio, the spectral exponents of the various primary components , and the 
secondaries from iron, that we can understand cosmic ray acceleration and 
propagation in a self consistent manner. Ws believe that with these measure- 
ments now being reported we can begin to separate acceleration and pro- 
pagation effects and that this represents an important step forward in 
cosmic ray research. 
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In this paper our results on high energy conposltlon are presented 
and their Implications are discussed. The preliminary analysis presented 
earlier (Ormes et al. , 1971, hereafter called OBR) was not sufficiently 

z' 

detailed or sophisticated to reveal all of these spectral differences. 

The results we present here are greatly improved by the use of a multidimen- 
sional charge analysis with more efficient background rejection, and a 
more accurate energy determination. Complex couplings between the charge, 
energy and trajectory information have been taken into account and will be 
discussed . 

II. EXPERIMENTAL DETAILS 

The balloon-borne instrument shown in figure 1, consisted of three 
major components, a charge measuring module, a spark chamber for determining 
particle trajectories, and an ionization spectrometer for measuring total 
energy. The charge of an incoming particle was determined using two plastic 
scintillators, a Lucite Cerenkov counter, and a Csl mosaic. Each detector 
had a sensitive area of 50 cm x 50 cm. The spark chamber was a digitized 
wire to wire chamber with four perpendicularly oriented planes, each with 200 
wires spanning 50 cm. For each particle four (x,y) measurements were 
available. The trajectory frcm the spark chamber was used to eliminate 
the dispersion in pulse height due to the variation of response over the area 
of the detector and due to the angle of incidence of the particles. A 
discussion of these corrections and details of the instrumentation have 
been presented previously (Ormes and Balasubrabmanyan, 1970.) 

The charge-module was followed by an electron cascade section consisting 
of a sandwich containing 12 tungsten plates and plastic scintillators. 
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Ihls section was designed to separate electrons from protons, and its 
thickness was 11 radiation lengths or about 1/2 a proton interaction mean 
free path (mfp). Results on the electron component are being published 
separately (Silverberg et al., 1973). The electron cascade section was 
followed by a nuclear cascade section consisting of 7 modules of iron, 
each 1/2 mfp thick. Each module had three plastic scintillators uniformly 
distributed inside the iron and viewed by a single photomultiplier tube at each 
of two opposite ends . Since a sample of the number of particles in the cascade 
was taken every 1.5 radiation lengths, the fluctuations due to low energy 
electron cascades were minimized. All the 23 detectors were pulse height 
analyzed and had dynamic ranges of 10 4 . 

The payload weighed 6000 lbs and it was flown on a mylar scrim 
balloon of 26 x 10^ ft 3 vdlume from Holloman Air Force Base in southern 
New Mexico on November 11, 1970. Useful data was obtained at a ceiling altitude 
of 7.4 gm/cm 2 for 14.4 hours. The geometric factor for the data presented 
was 725 cm^ s ter and the live time was 21.7 x 10^ sec. Particles at 
zenith angles greater than 25° were excluded, yielding a total exposure 
factor of 1420 cm2 ster sec. 

III. CHARGE DETERMINATION & BACKGROUND REJECTION 

The charge analysis is extremely critical because of its coupling with 
the energy determination and because of the background problems. As we shall 
see below, the response of the spectrometer depends upon energy/nucleon, 
and in order for this to be determined from the total observed energy 
deposit, the mass of the particles must be known. Furthermore, below 
about 2 GeV/nucleon, the Cerenkov (C) response is energy dependent, and so 
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the charge analysis depends upon the energy/nucleon. An interative pro- 
cedure is followed. First the charge is determined assuming that the C 
response is energy independent. From the charge a mass is assumed, and 
energy is determined. This energy is used to correct the predicted C 
response and the charge determination is rechecked. If it shifts, then 
the energy is redetermined. The scintillators at low energy show a slight 
increase due to the inverse square dependence of ionization loss on velocity, 
but since the effect is less than 5%, no correction was made. 

We consider first the details of the charge measurement. 

The charge was measured by the four detectors of the charge module. 

Ihe unique identification of a charge and the rejection of background 
depended upon two factors. The charge as measured by all four detectors 
had to be consistent within errors, and the trajectory of the particle 
had to be well defined and lie within the telescope geometry. When these 
conditions were met, pulse heights were corrected for geometrical variations 
in response using the trajectory data. 

Since the study included singly charged particles and iron nuclei as 

well as the ligat and medium nuclei, the tracks of the incident particles 

had to be determined over an ionization range of more than 600. With the 

o 

knock-on probability increasing as Z , heavier nuclei were invariably ac- 
companied by knock-on electrons which caused confusion in determining the 
track. A conputer algorithm was developed which detected the tracks of 
heavy (Z >3) nuclei with an efficiency estimated as 95%. 

Since the probability of a spark not forming along the track of the heavy 
nucleus was negligible, the main problem was to pick the correct trajectory 
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from a number of possible trajectories. The algorithm tried to fit straight 
trajectories through all combinations of sparks in the four decks of the 
wire grid spark chamber. The best fit to a straight line was selected with 
extra weighting given to sparks in which more than one adjacent wire partici 
pated. The mean number of sparks per event per plane increased from about 
2.5 at helium to about 5 at neon. This saturated the multiple sparking ef- 
ficiency of the chamber and so above neon there was no continuing increase 
in the number of sparks due to delta rays. Ey selecting tracks for the 
presence of sparks in which 2 or more wires participated, most of the pos- 
sible spurious tracks were eliminated. The trajectories determined by this 
algorithm were projected to their exit point in the iron spectrometer and 
were found highly correlated with sudden drops in the module signals. We 
conclude that "chamber inefficiency" is not a cause of missing a track, 
but the delta rays lead to seme uncertainty in track location. In order 
to check that this is not a problem, the zenith angle distributions of all 
nuclei have been checked to see that they correspond to expected distri- 
butions. If angular inaccuracies were more than a few degrees, the distri- 
butions would be broadened, and this was not observed. From this we esti- 
mate that o< 2° for the zenith angle measurements. 

Our confidence in the tracks is further enhanced because we have sub- 
sequently flown an experiment in which the number of decks was increased 
form 4 to 8. From an analysis of this new data it is clear that the tra- 
jectories picked by the 4 deck algorithm correspond closely to the 8 deck 
tracks . 



-7- 


TABLE 1 

Criteria for Selecting Simple Events: 

1. PARTICLE TRAJECTORY INTERSECTS COINCIDENCE SCINTILLATORS 

2. AIL 4 PLANES CONTAIN SPARKS 

3. LEAST SQUARES PIT OF DATA TO STRAIGHT LINE <1.5 WIRES 

4. THERE IS NO SECOND TRAJECTORY IN THE CHAMBER SATISFYING 1, 2 AND 3 
ABOVE WHICH HAS A DIFFERENT ZENITH ANGLE (> 5°) AND/OR POSITION (> 0.5in) 

Following the algorithm, events were divided into two classes, called 
simple and complex. Simple events had a single trajectory which passed all 
of the tests listed in table 1. The failure of any one would have classified it as 
complex. Only simple events were used in the analysis, and possible errors 
introduced by this procedure will be discussed below. 

The selection of simple events served to eliminate a large fraction 
of the background events. These background events were most abundant at 
Li and fell off rapidly with increasing charge. In particular, the spark 
chamber removed both atmospheric showers which could pass through the charge 
module unaltered and events caused by interactions in the spectrometer 
which triggered the experiment with back-scattered particles. The remaining 
background is eliminated by the charge determination procedure. 

To obtain the charge, the raw pulse heights were corrected for analyzer 
non- lin earities and for zenith angle and spatial variations in response. Due 
to the dispersions in the detectors, the observed pulse heights were dis- 
tributed around the centroids characteristic of each charge in the four 
dimensional representation. Single and two dimensional pulse height 
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distributions were constructed from corrected pulse heights, and the centroid 
and resolution were found for each charge and for each detector. A distance 
of each event from the charge centroid was determined by calculating the root 
mean square distance from the centroids in units of the resolution for the 
appropriate charge. For each particle, the charge was assigned from the 
centroid which gave the minimum distance. Particles were rejected if 
their average distance was more than about FWHM/2 away from the centroid. 

The correction for particles lost in this selection procedure will be dis- 
cussed below. 

In order to demonstrate the selectivity of the method, a one dimensio nal 
charge distribution was constructed. A charge was determined as the mini- 
mum in the parabola formed by the distance to the three closest centroids. 

(If for example, an oxygen nucleus passed through two detectors, lost a 

proton and became a nitrogen, its distance to the oxygen and nitrogen cen- 
troids would be minimized and its charge from the parabola fit would be 

approximately 7-5.) The charge histogram for charges 3 to 8 shown in 
figure 2 was constructed in this manner. This analysis has greatly im- 
proved the signal to noise ratio of the data over the earlier analysis 
(OBR) in which only two detectors were used for charge Identification. 

Individual charges are well resolved. The good peak to valley ratio 
shows that background rejection has been quite effective. In the iron 
range of the spectrum az = +1 charge unit, and so iron is separated 
from its fragmentation products. By far the most background exists in the 
Li and Be region, which is characteristic of experiments deep in the atmosphere, 
and so for these elements the spectra are likely to still have background 
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contamination . The problem Is that atmospheric showers from proton or 
helium interactions can produce 9 particles which in turn can pass through 
the charge module in a closely collimated bunch. Since the charge module 

contains only 0.16 radiation lengths of material, little cascade develop- 
ment takes place there and so the shower often passes through the charge 
module essentially unaltered. 

Let us consider the effects of the background rejection. In order 
to insure that background is minimized, strict criteria have been placed 
on the consistency of the pulse height measurements. This may mean that 
some events, especially those which interact near the top of the spectro- 
meter and produce backscatter,may be rejected. It is believed that this 
affects the acceptance of nuclei in the range of the light and medium 
nuclei, because the backscattered energy can be some 10 or 20 times mini- 
mum. However, this probably does not affect the heavier nuclei. It is 
difficult to identify these events in the data. Pulse heights are either 
totally inconsistent or the spark chamber trajectory is complex and a 
track cannot be found. Events which lie outside the error range but which 
have simple tracks could represent at most a 20$ correction to the data. 

The point at which this final cut (FWKM/2) is taken is such that the L/M 
ratio and the individual charge composition agrees with the published 
literature. Under these conditions, we do not believe our spectra for B 
and all higher charges are affected by background contamination. However, 
good events are undoubtedly rejected by our strict criteria. This problem 
will be the subject of further study, both with more detailed analysis and 
with an improved version of the spark chamber in which multiple trajectories 
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can be seen. 

To check that the spark chamber event classification scheme was not 
by itself biasing the data we analyzed all the complex events assuming 
they had trajectories at the most probable response angle of the telescope. 

To select these events they were required to deposit energy in the spectro- 
meter modules so that their zenith angles probably did not exceed 25°. These 
data were then subjected to the same charge analysis as used for the simple 
events. The data for all the classifications of conplex events listed in 
table 1 were consistent with being background. It is interesting to note 
that the number of "iron" group nuclei identified as complex is consistent 
with the expected number of interactions of iron nuclei in the charge 
module. For carbon nuclei, an upper limit of 18 % can be put on the losses 
to 'well identified complex events and to 25 % if marginally identified 
nuclei are allowed. The most probable correction is in the range of 5 to 
10 % . As discussed above, many carbon induced background events may be 
conplex, but they do not appear to have consistent pulse height measurements. 
We do not believe that this is due to delta rays in the chamber but due to 
interactions in the material of the spectrometer. 

IV. ENERGY MEASUREMENT 

The incident nuclei interacted in the ionization spectrometer and dis- 
sipated their energy via nuclear and electron cascades from tt° decay. The 
light signal seen by the photomultiplier tubes in the spectrometer was re- 
lated to the energy of the incident nucleus (Murzin, 1971, Akimov et al., 
1970, and Jones et al., 1969.) 

Most of the incident heavy nuclei interacted in the tungsten module. 
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and so a first approximation to the energy of the incident particles was 
obtained by taking the sum over all the modules of the number of equivalent 
muons in each module times the energy loss per module for relativistic muons. 
Using this energy, particles which pass through the entire spectrometer were 
grouped in bins and used to construct average integral cascade growth curves. 
These were done separately for the various charges and plotted on an energy/ 
nucleon basis . Two samples for carbon and iron are shown in figure 3 • 

Our earlier analysis (OBR) was based upon total energy and was not done 
separately for the different charges. Some differences between the curves 
can be seen at low energy and small depths. These differences are believed 
to be due to the importance of ionization loss at low energies. Above 2 
GeV/nucleon the curves are quite similar, and so all the data were combined 
on an energy/nucleon basis to get the series of growth curves shown in 
figure 4a. 

The curves could be expressed as: 

E = E (1 - exp(nA)) (1) 

no 

where Eh is the measured energy using n modules and E Q is the energy of 
the incident particle. ME q ), the number of modules necessary for absorbing 
(1-1/e) of the incident energy, was energy dependent. Ih figure 4b the 
energy dependence of X (E q ) is shown to increase from 1 to about 6 or 8 
modules over the energy range 1 to 100 GeV/nucleon. These empirically 
fitted curves also give the energy E c as a function of the measured energy. 
The correction used to obtain Eq from the data is shown in figure 4c. The 
spectral exponents derived using the two extreme values of the correction 
shown in the figure differ by less than half the statistical uncert ain ty. 
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The cascades developed by heavy cosmic ray nuclei in the spectrometer 
were compared with the Monte Carlo calculations of Jones (1971) . Our cas cad e 
curves agreed, and showed the energy deposit per incident nucleon was essen- 
tially indpendent of the charge of the nucleus. 

These curves also were used to find the energy of particles escaping 
out the sides of the spectrometer. These data can be used to double the 
statistical significance of our results. The data on nuclei exiting through 
the sides gave spectra 0.1 to 0.2 steeper than particles going through the 
bottom o However, the differences in spectral exponents between the dif- 
ferent nuclear species were approximately the same and the statistical 
significance was enhanced by including these particles. 

The final correction to the data was for the energy going into nuclear 
disintegration which is not observed by the spectrometer. This energy goes 
primarily into nuclear excitation, rest energy of the particles produced, 
and neutrons (Murzin, 1971) . This is perhaps the most uncertain portion 
of the analysis. We have compared the response of the spectrometer to the C 
response below 2 Ge V/i lucieon where the C varies sufficiently, and the two 
methods agree quite well. Most of these particles step in the spectrometer, 
and the correction made for nuclear disintegration energy is insignificant. 

At higher energy, we know from proton calibrations at 10 - 20 GeV/nucleon 
(Whiteside et al . , 1973) that 30 +5% of the energy was unobservable. This 
is in agreement with Monte Carlo calculations (Jones, 1971) which have now 
been compared with many calibration runs with different particles incl udin g 
C and O nuclei at 2.1 GV/c at the Bevatron at Berkeley. Using the Monte 
Carlo data we correct for a fraction of nuclear disintegration energy which 
decreases logarithmically from 50% at 1 GeV/nucleon to 20% at 100 GeV/ 
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nucleon. This represents a significant change from our earlier analysis 
(OBR) in which our corrections were based upon much less information and 
results in somewhat flatter spectra than those reported previously. 

Uncertainties in these corrections leave us with a possible systematic 
error in our spectral exponents of perhaps +0.1. Vfe believe, however, 
that the relative spectra are accurate to within the statistical uncertainties. 
The spectral exponents derived using only the energy deposited in the 
spectrometer are quite close to the spectra obtained after all corrections 
have been applied, giving us confidence that our corrections are not intro- 
ducing distortions in the shape of the energy spectra. 

Careful attention to the possible effects due to electronic non- 
linearities and detector response ncn-linearities convince us these effects 
make a negligible contribution. It is planned to check the response of the 
spectrometer against the response of a gas Cerenkov detector at energies above 
10 GeV/nucleon in the near future. This will provide us with an additional 
handle cn the nuclear disintegration correction. 

V. RESULTS 

The relative abundance observed at the balloon altitude has to be 
corrected for nuclear interactions in the matter in the telescope (4.5 
gm/cm 2 ) and for interactions in the atmosphere (7.4 gm/an 2 ) in order to 
get an accurate idea of the relative abundance distribution of the primary 
cosmic radiation. This requires a knowledge of all the fragmentation 
parameters and their energy dependence. Information on nuclear fragmentation 
parameters is still incomplete (Cleghom et al. , 1968, von Rosenvinge, 1970). 
Direct experimental determinations of fragmentation cross-sections using 
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protons incident on high Z targets have shewn that cross-sections are 
constant to +10% beyond 1 GeV/nucleon (Shapiro et al. , 1971). Though 
all fragmentation interactions have not been studied completely, vari- 
ation of cross-section with energy seems to be unlikely. 

In extrapolating the observed flux to the top of the atmosphere, 
we have used the concept of the absorption length Aj_ = / jjb (1-P^Nj/N^) 

where is the interaction mean free path, P^j are the fragmentation 
parameters for the production of the 1 th nucleus from heavier com- 
ponent, ISt and Nj are the abundances of i^ 1 and components respectively. 

We have used the Aj_ from Webber et al., (1972) to extrapolate the 
observed flux to the top of the atmosphere. 

Table 2 gives the integral flux of the different groups of nuclei at 
the top of the atmosphere. Also given are the results from other observers 
for comparison . 

TABLE 2 

Flux from Different Experiments in Particles/fri 2 -sec-ster 


Group of 
nuclei 

GSrC (>4.5 GV) 

Webber et al. 

Smith et al. 
(>5 GV) 

Mewaldt et al. 
(>5 GV) 

C+O 

3.45 + .7 

4.66 + 0.14 
(>4.35 GV) 

4.13 + .07 


10<Z<14 

1.08 + .2 


1.5 + .03 


15<Z<23 

.30 + .06 


.36 +.03 


20<Z<28 

.44 + .09 

.489 +.02 
for Z>17 
(>4.10 GV) 


.417 + .03 


Our carbon and oxygen fluxes are low compared to those of other workers, 
whereas there are no severe discrepancies for the other nuclei. We believe 
this is related to the background rejection we have discussed above. The 
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selection results in our losing seme genuine events up through oxygen nuclei. 
The correction for this loss is estimated to be at most 25% for carbon 
and oxygen and only a few percent for iron. This leaves us with a 35% 
discrepancy which is probably due to events interacting near the top of the 
spectrometer which lock like background in the charge module. 

The energy spectra of various nuclei and groups of nuclei are shown 
in figure 5. Below 2.0 GeV/nucleon the spectral shape is determined by 
the geomagnetic cutoff and was time dependent because the balloon drifted 
from 5.0 GV to 3.2 GV cutoff during the course of the flight. The data were 
fitted to power law spectra of the form 

dN/dE = k/E^ (E in GeV/nucleon) (2) 

above 2.0 GeV/nucleon. This energy was chosen in order to be free of 
geomagnetic cutoff effects. 

The data from all the nuclei are summarized in table 3. The expo- 
nents have been obtained using a least squares fitting technique. Five energy 
bins per decade were used giving a bin width comparable to the energy res- 
olution . The table also gives the 2 and the degrees of freedom for the 
fitted spectrum and the number of counts on which it is based. Vie believe 
that the difference in the spectral indices of the different species are 
known with greater confidence than the absolute value of the exponents of 
any nuclear species due to the systematic effects discussed. 

The H and He spectra are from data taken above 50 GeV total energy 
(Ryan et al. , 1972) . Seme of the differences in exponent between these 
and the heavier nuclei reflect a general flattening of the H and He spectra 
below 10 GeV/nucleon, probably due to solar modulation. 



TABLE 3 


Nuclei 

Y 

Z 

Degrees of 
Freedom 

Number of Events 
Included 

Number of Particles 
Above 20GeV/Nucleon 

Correction Factors 
For Flux* 

H 

He 

Li 

2.75+. 03 
2.77+.05 
2.28+. 15 

3.2 

6 

122 

2 

(1.25) (1.12) (1.13) 

Be 

2.6 +.2 

4.1 

4 

121 

0 

(1.25) (1.14) (1.13) 

B 

2.76+. 13 

2.2 

5 

229 

3 

(1.25) (1.16) (1.13) 

C 

2.52+. 06 

13.2 

7 

881 

25 

(1.25) (1.17) (1.27) 

N 

2.73+.11 

7.3 

6 

299 

5 

(1.25) (1.19) (1.24) 

1 

0 

2.57+. 06 

7.2 

7 

783 

25 

h* 

(1.25) (1.21) (1.30) T 

B+N 

2.77+.08 

4.1 

6 

528 

8 

(1.25) (1.17) (1.19+.05) 

C+O 

2.56+. 04 

14.0 

7 

1664 

50 

(1.25) (1.19) (1.27) 

10-14 

2.44+. 07 

.8 

7 

529 

18 

(1.15) (1.25) (1.32) 

15-23 

2.1 +.2 

1.3 

5 

140 

3 

(1.0) (1.33) (1.38) 

Fe group 

2.0 +.14 

2.4 

5 

99 

5 

(1.0) (1.45) (1.66) 

all Z>3 

2.55+. 03 

16.4 

8 

3207 

88 

(1.2) (1.2) (1.27) 


Geometric factor = 1420cm 2 ster sec 

*(correction for events outside identification limit) (correction for interactions in charge module) 
(correction for atmospheric attenuation) 
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The exponents of the energy spectra of the secondary nuclei B and 
N are larger than those of C and 0. The exponents of the heavy nuclei 
are considerably flatter. The spectra of Li and Be have been recovered 
from a region of very large background, and so our confidence in these 
spectra is much less. Their flatter spectra are characteristic of spectra 
with background contamination. 

The 2 values indicate that for the most part power laws are reason- 
able representations of the data. The spectral exponents agree well with 
those of Smith et al. , (1973) using a magnetic spectrometer. 

VI. DISCUSSION 

Our result that the B+N spectrum is steeper than the spectrum of C+0 is 
qualitatively consistent with the differences between L and M nuclei pre- 
viously reported: Webber and Ormes (1967) , smith et al., (1973) , and Webber 

et al., (1973). The data all indicate the spectra of secondary nuclei are 
steeper than those of primaries in agreement with Smith et al., (1973) 
and Juliusson et al., (1972) . We believe that the decrease in the secondary/ 
primary ratio is gradual and represents a power law difference of about 0.2 
in spectra, and does not indicate the precipitous drop seen by Juliusson 
et al. At these energies the fragmentation parameters are expected to be 
independent of energy, and so it is natural to follow Smith et al., in 
interpreting this result in terms of an energy dependence to the residence 
tine in the interstellar medium. In this case the lifetime probably decreases 
by a factor of 2 or 3 between 1 and 100 GeV/nucleon making the leakage 
length about 2 gnv/an 2 at 100 GeV/nucleon. If the leakage length does decrease 
faster than this, i.e. to a few tenths of a giVan 2 by 100 GeV/nucleon (Webber 
et al., 1973) then there will be severe difficulties with galactic models in 
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explaining the isotropy of underground muons from cosmic rays at 
10 12 eV (Elliot et al., 1970, Audouze and Cesarsky, 1973). 

WO believe (Ormas and Balasubrahmanyan, 1973, Ramaty et al. , 1973) 
that the difference between the Fe and C+O spectra is too large to be 
explained as a propagation effect as attempted by Webber et al., (1973) . 
Ibis difference in interpretation is in part related to experimental 
differences in the data from the secondaries produced by spallation of iron 
It can only be resolved by further experimentation . If the differences 
are related to the acceleratd.cn mechanism, it may be due to either separate 
sources (either location or type) or to Z dependent forces on the particles 
By measuring the spectra of Ne, and Si we should be able to resolve this 
question. Our spectra of the 10 to 14 group are consistent with having 
the same spectrum as the C+O within statistical uncertainties. 

The spectral exponents of secondaries are steeper than those of the 
primaries in these energy ranges, which suggests an energy dependent break- 
down of the trapping of these nuclei either in source regions or in the 
interstellar medium. If the flatter iron spectrum can be explained by 
propagation, the spectral exponent should steepen at high energies as 
the leakage length becomes shorter. 

The results reported in this paper represent the first use of an 
ionization spectrometer for the study of the energy distribution of heavy 
nuclei. Due to the limited exposure in the balloon flight, the data at 
high energy have statistical limitations. In the energy range over which 
most of the data has been collected (2 to 10 GeV/nucleon) there have been 
appreciable corrections due to the nuclear disintegration energy not seen 
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by the spectrometer. So far we have not been able to calibrate the re- 
sponse of the spectrometer to complex nuclei by another independent 
technique. we are planning to get seme checks on the energy measurement 
using a gas Cerenkov counter in a future flight. But in spite of the 
limitations of this data, it is already clear that there are many inter- 
esting changes in the composition of cosmic rays at high energies. These 
results point to the existence of a fertile field for experimental 
observations in the energy region beyond a few GeV/nucleon . For example, 
new results from the Chicago group (Juliusson and Meyer, 1973) indicate 
that the C/0 ratio also varies above 30 GeV/nucleon. 

Vfe should also remark that the other techniques which work in the range 
1 to 100 GeV/nucleon, such as gas Cerenkov detectors and magnetic spec- 
trometers have their limitations , and minor differences in the results are 
to be expected. In principle the spectraneter techniques should be best 
suited to the energy range 10 to 10 4 GeV/nucleon where particles are so rare 
that satellite observations of a year or more duration will be required to 
obtain statistically significant results. 

Future cosmic rays experiments in satellites should be able to resolve 
differences in interpretation of the data by being able to detect anisotropies, 
slight spectral discontinuities, and other small differences in carposition 
which may be missed by balloon experiments with their limited observing times. 
In any case we feel that the high energy carposition experiments in the several 
100 GeV region will give information crucial to the understanding of cosmic 
ray sources, acceleration phenomena, and propagation. 
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FIGURES 


Figure 1: 
Figure 2: 

Figure 3: 

Figure 4: 


Figure 5: 


Schematic diagram of the balloon flight experiment. 

Charge histogram obtained from 4 parameter charge analysis. 

Pulse height scale is normalized to 36 for carbon. 

Integral cascade growth curves for carbon and iron nuclei. 

Curves are similar on an energy/nucleon basis. 

Integral cascade growth curves for all nuclei are fitted 
to give the absorption length and thus the energy incident 
cn the spectrometer . 

Part a) Integral growth curves for all nuclei. The crosses 
represent the asymptotic energy and the numbers 
next to the crosses are A - \ The number of events 
on which these average curves are based are also given. 

Part b) The absorption mean free path as a function of energy. 

Part c) The correction factor giving the incident energy 

from the measured energy as a function of the measured 
energy. The correction used was the more extreme 
curve. Using the less extreme curve changed the spectral 
exponents by less than the 1 a statistical uncertainty. 
Energy spectra of various nuclei and groups of nuclei. The most 
obvious spectral differences are between C+0 and Fe. Data below 
2 GeV/nucleon are affected by geomagnetic cutoff. 
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